Complex periodic potentials with real band spectra 
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This paper demonstrates that complex "PT-symmetric periodic potentials possess real band 
spectra. However, there are significant qualitative differences in the band structure for these poten- 
tials when compared with conventional real periodic potentials. For example, while the potentials 
V{x) = isin 2JV+1 (x) (N = 0, 1, 2, . . .) have infinitely many gaps, at the band edges there are 
periodic wave functions but no antiperiodic wave functions. Numerical analysis and higher-order 
WKB techniques are used to establish these results. 



For a quantum mechanical model having a periodic potential the Schrodinger equation is 

-il)"(x)+V{x)il){x)=Eip{x), (1) 

where the potential V(x) is periodic with period P: 

V(x + P) = V(x). (2) 

In conventional treatments of Eq. (|l|) the periodic potential V(x) is assumed to be real. Imposing the condition 
' that the wave function ip(x) be bounded leads to a real spectrum consisting of continuous bands separated by gaps. 
^ ■ There is an infinite number of bands and gaps, except for the special family of so-called finite-gap potentials such as 
the Lame potentials ||. 

In this note we extend the conventional analysis to include the case of complex periodic potentials. We find that 
complex periodic potentials having VT symmetry exhibit real band spectra, despite the non-Hermitian character of 
the Schrodinger equation (|l|). (Here, V represents parity reflection and T represents time reversal.) Potentials having 
OO ' this symmetry satisfy 



[V(-x)]* = V(x). (3) 



Examples of such potentials are isinx, isin 3 (a;), and e lx . In addition to the property that these potentials have real 
I . spectra, the band structure displays several novel features that are strikingly different from the case of real periodic 
'"O ' potentials. 

The work reported here was motivated by recent investigations of non-Hermitian 'PT-symmetric Hamiltonian models 
having real discrete spectra. One such class of models is defined by the Hamiltonian B 



H = p 2 +x 2 (ix) £ (e>0). (4) 



Despite the lack of conventional Hermiticity, the spectrum of this Hamiltonian is real, positive, and discrete; each of 
the energy levels increases as a function of increasing e. It has been observed that the reality of the spectrum is a 
consequence of VT symmetry, which is a weaker condition than Hermiticity. This observation has also been used to 
construct new classes of quasi-exactly solvable quantum theories [0 and to study new kinds of symmetry breaking in 
quantum field theory There have been many other instances of non-Hermitian VT- invariant Hamiltonians in 

physics. Energies of solitons on a complex Toda lattice have been found to be real |^. Hamiltonians rendered non- 
Hermitian by an imaginary external field have been used to study population biology || and to study derealization 
transitions such as vortex flux-line depinning in type-II superconductors 113). In these last two cases, initially real 
eigenvalues bifurcate into the complex plane due to the increasing external field, indicating the growth of populations 
or the unbinding of vortices. 

We begin by summarizing the standard Floquet analysis of the Schrodinger equation (|]J) for the case where V(x) is 
real and periodic We define a fundamental pair of linearly independent solutions u\{x) and U2(x) satisfying the 
initial conditions 

ui(0) - 1, ui(Q) - 0; 

u a (0) = 0, u' 2 (0) = l. (5) 

Any solution ip(x) to Eq. (Q) is a linear combination of u\{x) and U2(x). It is then a straightforward algebraic exercise 
to show that ip( x ) is bounded provided that the discriminant A(E), which is defined by 
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A(E) = Ul (P)+u' 2 (P), 



(6) 



satisfies the constraint 



2 < A(E) < 2. 



(7) 



To illustrate the features of the discriminant we consider a typical periodic potential, V(x) = sin(a;), for which the 
period P = 2tt. In Fig. [H we plot A(E) as a function of E. Note that A(E) is oscillatory and is well approximated 
by the function 2cos(7rvi5) for large E. The crucial feature of A(E), which cannot be seen from this plot, is that 
its graph crosses the lines ±2 an infinite number of times; each of the maxima of A(E) lies above 2 and each of the 
minima lies below —2. The regions of energy for which |A(2?)| < 2 are called bands and the regions of energy for 
which |A(i?)| > 2 are called gaps. The gap size decreases exponentially as a function of E. The band edges at which 
A(E) — 2 correspond to periodic solutions to Eq. (|l|), ip(x) = ip(x + P), and the band edges at which A(E) = —2 
correspond to antiperiodic solutions %j)(x) = —ip(x + P). 
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FIG. 1. The discriminant A(E) plotted as a function of E for the real periodic potential V(x) = sin(a;). Although it cannot 
be seen in the figure, all local maxima lie above the line A — 2 and all local minima lie below the line A = —2. The regions of 
energy E for which |A| < 2 correspond to bands and the regions where |A| > 2 correspond to gaps. There are infinitely many 
gaps and these gaps become exponentially narrow as E increases. 



Now consider the calculation of the discriminant for the case of a complex periodic potential V(x). In general, a 
complex periodic potential will have no bounded solutions because the discriminant is typically complex. However, 
for complex T-'T-symmetric periodic potentials, one can easily show that the discriminant A(E) is real when E is 
real. The VT symmetry is crucial here; for a potential that is not VT symmetric [one that does not satisfy Eq. (0)], 
the discriminant is complex for all values of E. Having established that complex 'PT-symmetric periodic potentials 
have real discriminants, we can then apply the criterion in Eq. (Q) to locate the real energy bands within which the 
corresponding wave function ip(x) is a bounded function. 

We have computed the discriminants for the class of complex 'PT-symmetric periodic potentials 



V(x) 



2N+1 



(x) (N = 0, 1, 2, ...). 



(8) 



In Figs. we plot the discriminants for the cases N = 0, 1, . . . , 5. While these plots superficially resemble Fig. [l] 
for large E, they exhibit new and intriguing features that are significantly different from the case of a real periodic 
potential. The most obvious new feature is the appearance for N > 1 of a local minimum of A(E) between —2 and 
2. Such a dip is rigorously forbidden in the case of real periodic potentials 

The most dramatic differences between the discriminants for the complex PT-symmetric periodic potentials in 
Eq. (||) and real periodic potentials cannot be easily seen in the figures. We have performed a careful numerical study 
of the local minima and maxima of A(E). Our study reveals that none of the local minima lies below —2. This shows 
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that there are no antiperiodic solutions if)(x) to the Schrodinger equation (|l|). Nevertheless, all of the local maxima of 
A(_E) lie above 2. Hence, there are an infinite number of band gaps in the spectrum and the band-edge wave functions 
are periodic. 
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FIG. 2. The discriminant A(E) plotted as a function of E for the complex periodic potential V(x) = isin(a;). On the scale 
of this figure for E > 2 it is not possible to see any difference between this figure and Fig. 1. However, although it cannot be 
seen in this figure, all local maxima lie above the line A = 2 and all local minima lie above the line A = —2. This behavior is 
distinctly different from the generic behavior in Fig. 1 exhibited by the real periodic potential V(x) = sin a;. 
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FIG. 3. The discriminant A(E) for the complex periodic potential V(x) = isin 3 (x). All local maxima lie above the line 
A = 2 and all local minima lie above the line A = —2. Note the appearance of a new local minimum near E = 0.5. 
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FIG. 4. The discriminant A(E) for the complex periodic potential V(x) = isin 5 (:r). All local maxima lie above the line 
A = 2 and all local minima lie above the line A = —2. 
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FIG. 5. The discriminant A(E) for the complex periodic potential V(x) 
A = 2 and all local minima lie above the line A = —2. 



isin 7 (x). All local maxima lie above the line 
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FIG. 6. The discriminant A(E) for the complex periodic potential V(x) = isin 9 (x). All local maxima lie above the line 
A = 2 and all local minima lie above the line A = —2. 
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FIG. 7. The discriminant A(E) for the complex periodic potential V(x) = isin (a;). All local maxima lie above the line 
A = 2 and all local minima lie above the line A = —2. 



To perform this numerical analysis it is necessary to locate the positions of the local minima and maxima of A(E) 
to extremely high accuracy. This can be done using WKB methods @. We take the energy E to be large (E » 1) 
and define a small parameter e by 

(9) 



E 

Then we make an exponential ansatz for the wave function ip( x ) 
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ip(x) — exp 



(10) 



n=0 



Substituting ip{x) m Eq. (|Io| ) into the Schrodinger equation (Q) gives a recursion relation for the functions Q n (x): 

l-[Q' (x)} 2 = 0, 
iQZ(x)-2Q' (x)Q[(x) = 0, 
iQ'Kx) - 2Q' (x)Q' 2 (x) - lQ[(x)] 2 - V(x) = 0, 



iQn-M -J2Q'A x )Q'n-j(x) = (n > 3). 



(11) 



The solution to these equations is 



Qo(x 
Qi(x 

Q2{x 
Qz{x 
Qa{x 
Q 5 (x 
Q 6 (x 



±x, 
0, 



i r x 
= T 2 1 dtv{t) > 

= -i{V(x)-V(Q)}, 
= ±l(v\x)-V'(0)- I dtV 2 (t)\, 



— [V"(x) - V"(0) - 2V 2 (x) + 2V 2 (0)} , 

T ~k ( V "'W ~ V "'^ ~ 5V ( X ) V '( X ) + 5V(0)V"(0) + I" dt [2V 3 {t) - V{t)V"(t)] 



and so on. Note that Q n {x] is normalized so that Q n (Q) = 0. In general, the formula for Q' n (x) is 



QU X ) = ±2 



3=1 



(n > 3). 



(12) 



(13) 



In order to obtain a WKB formula for the discriminant A(E) in Eq. we need to evaluate Q n (x) at x = P. The 
periodicity of the potential V(x) simplifies the results considerably; when n is odd, Q n (P) — and when n is even, 
only the integrals in Eq. ( |l2|) remain: 

Qo(P) = ±P, 



1 f P 

Q 2 (P)^T^ j dtV(t), 

Qa{P) = t\J^ dtV 2 {t), 

Qe(P) = ^ dt [2V\t) - V(t)V"(t)}, 



(14) 



and so on. 

The WKB formula for the discriminant is particularly simple when the potential V(x) is VT symmetric: 



A(E) = 2 cos 



1 OO 

-^e 2 "Q 2 »(P) 



n=0 



(15) 



One obtains the same formula for potentials that are real and symmetric under parity V . 

For the complex 'PT-symmetric potentials V(x) in (|J) the WKB formula for the discriminant in Eq. ([l5|) is 



A(E) = 2 cos 



2tt r(27V + 3/2)V^ 3 (2iV+ 1)T(2N + 1/2)^ 5 
T + 4(2iV+l)! 6 + 32(2iV)! 6 



(16) 
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A similar WKB formula exists for the real odd-parity potentials V(x) = sm 2N+1 (x): 

"2tt T(2N + 3/2)0? £ 3 (2JV + l)r(2iV + 1/2) y/n ^ 



A(E) = 2 cos 



4(2AT+1)! 32(2iV)! 



(17) 



We can illustrate the extreme accuracy of these WKB approximations by comparing them with numerical com- 
putations of the discriminant. For example, for V(x) = isin 3 (a;) at e = 0.2 (which corresponds to E — 25) we find 
numerically that A(25) = 1.9999960002. The first three orders of the WKB approximation taken from Eq. ( fl6|) give 
A(25) = 2.0, A(25) = 1.9999961447, and A(25) = 1.9999960046. Similarly, for V(x) = sm 3 {x) at e = 0.2 we find 
numerically that A(25) = 1.9999959937. The first three orders of the WKB approximation taken from Eq. ( fl7|) give 
the same values: A(25) = 2.0, A(25) = 1.9999961447, A (25) = 1.9999960046. 

Despite this impressive precision, the WKB formulas jlq ) and (|l7|) cannot be used directly to answer the crucial 
question of whether there are band gaps because these approximations to the discriminant A(E) never cross the values 
±2. The reason for this inadequacy of the WKB approximation is that the differences |A max — 2| and |A m i n + 2| are 
exponentially small when E » 1. Therefore, these differences are subdominant with respect to the WKB asymptotic 
series and are not accessible to any order in powers of e. Indeed, the WKB series can only provide information about 
quantities with an algebraically small error, and not an exponentially small error. We emphasize that the WKB 
approximation has this shortcoming only at the maxima and minima of the approximation. At other points any 
exponential discrepancy is completely negligible compared with algebraic errors. 

The WKB series is still an extremely useful ingredient in the numerical search for zeros of A(E) ± 2. (These zeros 
are the dividing points between bands and gaps.) Our procedure is first to find the energies at which there are maxima 
and minima of the WKB approximation to the discriminant and then to evaluate, with high numerical precision, the 
actual value of the discriminant in a tiny neighborhood of each of these points. By doing this we can determine 
whether or not the discriminant A(E) crosses the lines ±2. 

For the real periodic potentials V(x) — sin 2N+1 (x) our procedure confirms the rigorous theoretical result that every 
maximum of the discriminant lies above 2 and every minimum lies below —2. Consider, for example, the potential 
V(x) — sin(x). From Fig. 1, it is clear that the first maximum lies above 2. The first minimum occurs at E = 2.3138, 
where the discriminant has the value —2.0038787. The second maximum occurs at E = 4.0336, where the discriminant 
is 2.000007. Similar behavior is found for the other potentials in the class V{x) = sm 2N+1 (x). In stark contrast, for 
the potentials V(x) = i sm 2N+1 (x), while the maxima of the discriminant lie above +2, the minima of the discriminant 
lie above —2. Thus, for these potentials there are no antiperiodic wave functions. As an example, lengthy and delicate 
numerical analysis verifies that for the potential V(x) — ism(x) the first three maxima of the discriminant A(E) are 
located at E = 3.9664284, E = 8.9857320, and E = 15.99206621346. The value of the discriminant A(E) at these 
energies is 2.000007, 2.00000000000069, and 2.000000000000000000015. The first two minima of the discriminant are 
located at E = 2.1916 and E = 6.229223 and at these energies A(E) has the values -1.9953386 and -1.99999999527. 
Similar behavior is found for the other potentials in the class (Js|) . 

We conclude by pointing out that from the expressions for Q n {P) in Eq. (|IJ) the WKB series truncates if the 
potential is a polynomial in e lx . For example, for the complex 'PT-symmetric periodic potential V(x) — e 11 the 
WKB series in Eq. ( plf ) truncates after the first term because Q n {2ir) vanishes for n > 1. For this case the WKB 
approximation is exact and the discriminant is given by 



A(E) = 2cos(27rV£0. (18) 

One can verify this result directly by solving the Schrodinger equation (Q) for this potential exactly; the solution is a 
Bcssel function: ijj(x) = J 2 ^{2e lx / 2 ). 
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